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The rapid experimental progress of ultra-cold dipolar fermions opens up a whole new oppor¬ 
tunity to investigate novel many-body physics of fermions. In this article, we review theoretical 
studies of the Fermi liquid theory and Cooper pairing instabilities of both electric and magnetic 
dipolar fermionic systems from the perspective of unconventional symmetries. When the electric 
dipole moments are aligned by the external electric field, their interactions exhibit the explicit 
d r 2_ 3z 2 anisotropy. The Fermi liquid properties, including the single-particle spectra, thermody¬ 
namic susceptibilities, and collective excitations, are all affected by this anisotropy. The electric 
dipolar interaction provides a mechanism for the unconventional spin triplet Cooper pairing, which 
is different from the usual spin-fluctuation mechanism in solids and the superfluid 3 He. Further¬ 
more, the competition between pairing instabilities in the singlet and triplet channels gives rise to 
a novel time-reversal symmetry breaking superfluid state. Unlike electric dipole moments which are 
induced by electric fields and unquantized, magnetic dipole moments are intrinsic proportional to 
the hyperfine-spin operators with a Lande factor. Its effects even manifest in unpolarized systems 
exhibiting an isotropic but spin-orbit coupled nature. The resultant spin-orbit coupled Fermi liquid 
theory supports a collective sound mode exhibiting a topologically non-trivial spin distribution over 
the Fermi surface. It also leads to a novel p -wave spin triplet Cooper pairing state whose spin 
and orbital angular momentum are entangled to the total angular momentum J = 1 dubbed the 
J-triplet pairing. This J-triplet pairing phase is different from both the spin-orbit coupled 3 He-I? 
phase with J = 0 and the spin-orbit decoupled 3 He-T phase. 

PACS numbers: Keywords: electric and magnetic dipolar interactions, anisotropic Fermi liquid theory, SO 
coupled Fermi liquid theory, p-wave triplet Cooper pairing, time-reversal symmetry breaking 


I. INTRODUCTION 

Dipolar interactions have become a major research 
focus of ultra-cold atomic and molecular physics. For 
bosonic atoms with large magnetic dipolar moments (e.g. 
52 Cr), their magnetic moments are aligned in the Bose- 
Einstein condensates in which the anisotropy of the dipo¬ 
lar interaction is manifested^—. On the other hand, 
the synthesis and cooling of both fermions with elec¬ 
tric and magnetic dipolar moments give rises to an even 
more exciting opportunity to explore novel many-body 
physics^c— . The quantum degeneracy of the fermionic 
dipolar molecules of 40 K 87 Rb has been realized^—. 
These molecules have been loaded into optical lattices in 
which the loss rate is significantly suppressedii^. The 
chemically stable dipolar molecules of 23 Na 40 K have been 
cooled down to nearly quantum degeneracy with a life¬ 
time reported as 100ms near the Feshbach resonance^. 
The quantum degeneracy of fermionic atoms with large 
magnetic dipole moments has also been achieved for the 
systems of 161 Dy with — — and 167 Er with 7 

, which are characterized by the magnetic dipolar 
interaction. 

Electric and magnetic dipolar fermions exhibit novel 
many-body physics that is not well-studied in usual 
solids. One of the most prominent features of the elec¬ 
tric dipolar interaction is spatial anisotropy, which is 
markedly different from the isotropic Coulomb interac¬ 
tion in solids. In contrast, the magnetic dipolar interac¬ 


tion remains isotropic in unpolarized systems. More im¬ 
portantly, it exhibits the spin-orbit (SO) coupled feature, 
i.e., the magnetic dipolar interaction is invariant only un¬ 
der the simultaneous rotation of both the orientations of 
magnetic moments and their relative displacement vec¬ 
tors. These features bring interesting consequences to 
the many-body physics of dipolar fermions. 

Rigorously speaking, so far there are still no permanent 
electric dipole moments having been discovered yet at 
the levels of the elementary particle, atom, and molecule. 
For example, for a hetero-nuclear dipolar molecule, even 
though at an instantaneous moment, it exhibits a dipole 
moment, while it averages to zero in the molecular rota¬ 
tional eigenstates. External electric fields are needed to 
polarize electric dipole moments, which mixes rotational 
eigenstates with opposite parities. However, the dipole 
moment of these mixed states is unquantized, and, thus 
the electric dipole moment is a classic vector. When two 
dipole moments are aligned, say, along the 2 -axis, the 
interaction between them is spatially anisotropic, which 
not only depends on the distance between two dipoles, 
but also the direction of the relative displacement vector. 
Nevertheless, this anisotropy exhibits an elegant form of 
the spherical harmonics of the second Legendre polyno¬ 
mial, i.e., the d r 2 _ 3z 2 -type anisotropy22r— . This elegant 
anisotropy greatly simplifies the theoretical study of the 
novel many-body physics with the electric dipolar inter¬ 
action. 

The electric dipolar interaction results in an 
anisotropic Fermi liquid state, which exhibits differ- 
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ent single-particle and collective properties from those 
of the standard isotropic Fermi liquid theoryZir— . 
The shape of the Fermi surface exhibits anisotropic 
distortion s 31 ' 32 ! 34 ' 40 ! 43 . The anisotropic dipolar interac¬ 
tion mixes different partial-wave channels, and thus the 
usual Landau interaction parameters in the isotropic case 
should be generalized into the Landau interaction matrix 
with a tri-diagonal structure, which renormalizes ther¬ 
modynamic susceptibilitie s 33 ' 34 . The dispersion of the 
collective zero sound mode is also anisotropic: the zero 
sound mode can only propagate in a certain range of the 
solid angle direction, and its sound velocity is maximal 
if the propagation direction is along the north or south 
pole s 34 ' 36 . 

The anisotropy of the electric dipolar interac¬ 
tion also results in unconventional Cooper pairing 
symmetrie s 28 ' 29 ' 44- —. The electric dipolar interaction 
is neither purely attractive nor purely repulsive. The 
partial-wave analysis shows that the most attractive pair¬ 
ing channel is p z -like, which naturally gives rise to a new 
mechanism to unconventional pairing symmetry. Conse¬ 
quently, for the single component case, the pairing sym¬ 
metry is mostly of p z -like slightly hybridized with even 
higher odd partial wave component s 28 ' 29 ' 44 ' 45 . The pair¬ 
ing structure of the two-component dipolar fermions is 
even more interesting, which allows both the s + d -wave 
channel singlet and the p z -wave triplet pairing s 50 ' 52 ' 53 ' 55 . 
The dipolar interaction induced triplet pairing is to first 
order in interaction strength. In comparison, the spin 
fluctuation mechanism in solid state systems (e.g. 3 He 
and Sr 2 Ru 04 ) is a higher order effect of interaction s 56 ' 57 . 
The singlet and triplet pairing symmetries can coexist in 
two-component electric dipolar fermion systems. Only 
when their relative phase angle is ±lj, the resultant pair¬ 
ing is unitary— This gives rise to a novel and very 
general mechanism to a spontaneous time-reversal (TR) 
symmetry breaking pairing state. 

Next we discuss the novel feature of the magnetic dipo¬ 
lar fermion s 43 ' 58- —. The magnetic dipolar interaction 
is very complicated to handle in classic systems, which 
leads to a variety of rich patterns in real space. In com¬ 
parison, for the quantum degenerate Fermi systems, the 
existence of Fermi surfaces constraints the low energy 
degrees of freedom only around the Fermi surface. This 
feature greatly simplifies the theoretical analysis, and the 
exotic physics with non-trivial spin texture patterns lies 
in momentum space instead of real space. 

Typically speaking, the interaction energy scale of 
magnetic dipolar fermions is much smaller than that 
of the electric dipolar case. Nevertheless, conceptually 
they are still very interesting. Unlike the electric dipo¬ 
lar moment, the magnetic moment is proportional to 
the hyperfine-spin with the Lande factor, and thus its 
components are non-connnutative quantum-mechanical 
operator s 60 ' 62 . Magnetic dipole moments are permanent 
in the sense that they do not need to be induced by ex¬ 
ternal magnetic fields. In the absence of external fields, 
the unpolarized magnetic dipolar systems are in fact 


isotropic. Neither spin nor orbital angular momentum 
is conserved; nevertheless, the total angular momentum 
remains conserved by the dipolar interaction. Thus the 
magnetic dipolar interaction exhibits the essential feature 
of the SO coupling physics. Very recently, using electric 
dipolar moments to generate effective SO coupled inter¬ 
actions similar to that in the magnetic dipolar systems 
is also proposed in Ref. E3 by properly coupling mi¬ 
crowaves to molecular rotation eigenstates. 

The ordinary SO coupling in solids is a single-particle 
effect originating from the relativistic physics. In con¬ 
trast, in magnetic dipolar fermion system s 59 ' 61 ' 62 , the 
Fermi surfaces remain spherical without splitting in 
the absence of the external magnetic fields. Never¬ 
theless, this SO coupling appears at the interaction 
level, including the SO coupled Fermi surface Pomer- 
anchuk instabilitie s 59 ' 61 1 62 , and topological zero-sound 
wave modes exhibiting an oscillating spin distribu¬ 
tion of the hedgehog-type configuration over the Fermi 
surface— 

The magnetic dipolar interaction also induces novel 
Cooper pairing structures exhibiting the SO coupled 
natur o 60 ' 66 . Even in the simplest case of F = the 
magnetic dipolar interaction provides a novel and robust 
mechanism for the p -wave (L = 1) spin triplet (S = 1) 
Cooper pairing which arises from the attractive channel 
of the magnetic dipolar interaction. It turns out that its 
pairing symmetry structure is markedly different from 
that in the celebrated p -wave pairing system of 3 He: the 
orbital angular momenta L and spin S of Cooper pairs 
are entangled into the channel of the total angular mo¬ 
mentum J = 1, dubbed as the J-triplet pairing. In com¬ 
parison, the 3 He-I? phase is isotropic in which J = 0, 
while the A-phase is anisotropic in which J is not well- 
defined^. 

In this article, we review the recent progress of the 
novel many-body physics with dipolar fermions, such as 
the Fermi liquid properties and Cooper pairing struc¬ 
tures, focusing on unconventional symmetries. In Sect. 
El we review the anisotropy of the electric dipolar in¬ 
teraction, and the SO structure of the magnetic dipo¬ 
lar interactions, respectively, from the viewpoint of their 
Fourier components. In Sect. ED the anisotropic Fermi 
liquid theory of the electric dipolar fermions is reviewed. 
And the SO coupled Fermi liquid theory of the magnetic 
dipolar fermion systems is reviewed in Sect. IIVI The p z - 
wave Cooper pairing in the single and two-component 
electric dipolar systems and the TR reversal symmetry 
breaking effect are reviewed in Sect. |V] The SO coupled 
Cooper pairing with the J-triplet structure in the mag¬ 
netic dipolar fermion systems is reviewed in Sect. IVII 
Conclusions and outlooks are presented in Sect. I VIII 

Due to limit of space and also the view point 
from the unconventional symmetry, we are not able 
to cover many important research directions of dipo¬ 
lar atoms and molecules in this review. For example, 
the progress on topics of strong correlation physics with 
dipolar fermions^Zr— , the Feshbach resonance with dipo- 
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lar fermion s 52 i 53 , the synthetic gauge field with dipolar 
fermion s 70 ’ 71 , and the engineering of exotic and topologi¬ 
cal many-body states^- 4 . Some of these progresses have 
been excellently reviewed in Ref. [l7ll30j . The properties 
of dipolar boson condensations are not covered here ei¬ 
ther, and there are already many important reviews on 
this topi&b&^ i 30 i 75 i 76 . 


II. FOURIER TRANSFORM OF DIPOLAR 
INTERACTIONS 

In this section, we review the Fourier transformations 
of both the electric and magnetic dipolar interactions in 
Sect. Ill Al and Sect. IIIBI respectively. The anisotropy 
of the electric dipolar interaction and the SO coupled 
feature of the magnetic dipolar interaction also mani¬ 
fest in their momentum space structure. These Fourier 
transforms are important for later analysis of many-body 
physics. 


A. Electric dipolar interaction 


Without loss of generality, we assume that all the elec¬ 
tric dipoles are aligned by the external electric field E 
along the ^-direction, then the dipolar interaction be¬ 
tween two dipole moments i s- - > — 

2d 2 

V d (ri 2 ) = -t- -P 2 (cos0 12 ), (1) 

R 12 I 3 

where d is the magnitude of the electric dipole moment; 
fi 2 = r*i — f 2 is the displacement vector between two 
dipoles; di 2 is the polar angle of ri 2 ; P 2 { cos$ 12 ) is the 
standard second Legendre polynomial as 

P 2 (cos $ 12 ) = i(3 cos 2 012 - 1). 


The zeros of the second Legendre polynomial lie around 
the latitudes of 9q and ir — 9q with 


9 0 = cos 


1 

71 


55°. 


( 2 ) 


Within 6q < 9\ 2 < 7r — $o, the dipolar interaction is 
repulsive, and otherwise, it is attractive. The spatial 
average of the dipolar interaction in 3D is zero. 

For later convenience, we introduce the Fourier trans¬ 
form of the dipolar interaction Eq. [T] 


V d (q) = I d 3 re~^V d (r). 


( 3 ) 


same symmetry factor of the spherical harmonics. Sec¬ 
ond, since V d (r) decays with a cubic law, V d (q) should be 
dimensionless. 

If V d (fi 2 ) were isotropic, V d (q) would logarithmically 
depends on q. However, a more detailed calculation 
shows that actually it does not depend on the magni¬ 
tude of q. Let us introduce a short distance cutoff e that 
the dipolar interaction Eq. [l]is only valid for r > e, and 
a long distance cutoff R as the radius of the system. A 
detailed calculation shows that^i 

V d {q) = 8tt d 2 1 p 2 ( CO s6>g-) (4) 

{ qe qH ) 


where ji(x) is the first order spherical Bessel function 
with the asymptotic behavior as 



-1 

in 


as x —y 0; 
as x —> 00 . 


( 5 ) 


After taking the limits of qe —> 0 and qR -A + 00 , we 
arrive at 


V d (q) 


8ird 2 

3 


Pi (cos 9^). 


( 6 ) 


At q = 0, Vd{q = 0) is defined as 0 based on the fact 
that the angular average of the 3D dipolar interaction 
vanishes, thus V d is singular as q —> 0. Even in the case 
that R is large but finite, the smallest nonzero value of qR 
is at the order of 0(1). Thus, V d (q) remains non-analytic 
as q -A 0. 

An interesting feature of the above Fourier transform 
Eq. ©is that the anisotropy in momentum space is oppo¬ 
site to that in real space: it is most negative when q lies 
in the equatorial plane, and most positive when q points 
to the north and south poles. An intuitive picture is ex¬ 
plained in Fig. [TJ Consider a spatial distribution of the 
dipole density p(r), then the classic interaction energy is 


J J dridr 2 p(fi)p(r 2 )V d (fi - r 2 ) = 7- ^ \p{q)\ 2 

q 

X v d (q), (7) 

where Vo is the system volume. If the wave vector q is 
along the z-axis, then the dipole density oscillates along 
the dipole orientation, thus the interaction energy is re¬ 
pulsive. On the other hand, if q lies in the equatorial 
plane, the dipole density oscillates perpendicular to the 
dipole orientation, and thus the interaction energy is at¬ 
tractive. 


B. Magnetic Dipolar Interaction 


A lot of information can be obtained solely based on sym¬ 
metry and dimensional analysis. First, e lq ' r is invariant 
under spatial rotations, thus V d (q) transforms the same 
as V d (f) under spatial rotations. It should exhibit the 


Now let us consider the magnetic dipolar 
interaction^-^* 26 - 27 . Different from the electric dipole 
moment, the magnetic one originates from contributions 
of several different angular momentum operators. The 
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FIG. 1: The Fourier components of the dipolar interaction 
Vd(k). The left-hand-side is for k || z, and the right-hand-side 
is for k -L z. 


total magnetic moment is not conserved, and thus its 
component perpendicular to the total spin averages to 
zero. For the low energy physics below the coupling 
energy among different angular momenta, the magnetic 
moment can be approximated as just the component 
parallel to the spin direction, and thus the effective 
magnetic moment is proportional to the hyperfine spin 
operator up to a Lande factor, and thus is a quantum 
mechanical operator. Due to the large difference of 
energy scales between the fine and hyperfine structure 
couplings, the effective atomic magnetic moment below 
the hyperfine energy scale can be calculated through the 
following two steps. The first step is the Lande factor for 
the electron magnetic moment respect to total angular 
momentum of electron defined as jl e = <7j/.4bJ, where 
[jLb is the Bohr magneton; J = L + S is the sum of 
electron orbital angular momentum L and spin S; and 
the value of gj is determined as 

9 L+g s , g L ~ g s L(L + \)-S(S+1) 

9j - —o- 1 -o- rn i i\ -■ W 


Further considering the hyperfine coupling, the total 
magnetic momentum is defined jl = [Ib{, 9 jJ + 9il ) = 
g F F where gi is proportional to the nuclear gyromag- 
netic ratio and is thus tiny, and F is the hyperfine spin. 
The Lande factor g F can be similarly calculated as 


9f 


gj + gi gj - gi J{J + 1 ) - i{i + 1 ) 

2 2 F(F +1) 

g j( J(J + 1) — I (I + 1) \ 

2 V F(F + 1) /' 


(9) 


The magnetic dipolar interaction between two spin-F 
atoms located at fq and r-i is 


Fa/3;/3' a' (g) 


9f9 2 b 

y*3 



3(F aa i ■ r)(Fpp> ■ f) , 

( 10 ) 


where r = rq — r% and r = r/r is the unit vector along 
r. Similarly to the case of the electric dipolar interac¬ 
tion, the Fourier transform of Eq. [10] possesses the same 
symmetry structure as that in real spac o 43 ’ 59 


VapiP'a' (g) 


4tt g 2 F g 2 B 
3 


3 (F aa > ■ q){Fpp- ■ q) 



( 11 ) 


Again, it only depends on the direction of the momentum 
transfer but not on its magnitude, and it is also singular 
as q —> 0. If q is exactly zero, V a p-,p' a '(q = 0) = 0. 

In the current experiment systems of magnetic dipolar 
atoms, the atomic spin is very large. For example, for 
161 Dy, its atomic spin reaches F = and thus an accu¬ 
rate theoretical description of many-body physics of the 
magnetic dipolar interactions of such a large spin system 
would be quite challengin g -!*—. Nevertheless, as a theo¬ 
retical starting point, we can use the case of F = | as a 
prototype model which exhibits nearly all the qualitative 
features of the magnetic dipolar interaction s 43 ’ 60 . 


III. ANISOTROPIC FERMI LIQUID THEORY 
OF ELECTRIC DIPOLAR FERMIONS 

In this section, we will review the new ingredients of 
the Fermi liquid theory brought by the anisotropic elec¬ 
tric dipolar interactional"2£ i 41 i 42 ; including the single¬ 
particle properties such as Fermi surface distortions, and 
two-body properties including thermodynamic properties 
and collective modes. 

A general overview of the Landau-Fermi liquid the¬ 
ory is presented in Sect. IIII Al In Sect. IIIIB1 we re¬ 
view the dipolar interaction induced Fermi surface dis¬ 
tortions. The Landau interaction matrix is presented in 
Sect. IIII Cl and its renormalization on thermodynamic 
properties including Pomeranchuk instabilities are review 
in Sect. mm The anisotropic collective excitations are 
reviewed in Sect. IIII El 


A. A quick overview of the Fermi liquid theory 

One of the most important paradigms of the in¬ 
teracting fermion systems is the Landau Fermi liquid 
theory!!" 79 . The Fermi liquid ground state can be viewed 
as an adiabatic evolution from the non-interacting Fermi 
gas by gradually turning on interactions. Although the 
ideal Fermi distribution function could be significantly 
distorted, its discontinuity remains which still defines a 
Fermi surface enclosing a volume in momentum space 
proportional to the total fermion number, as stated by 
the Luttinger theorem. Nevertheless the shape of the 
Fermi surface can be modified by interactions. The low 
energy excitations become the long-lived quasi-particles 
around the Fermi surface, whose life-time is inversely pro¬ 
portional to the square of its energy due to the limited 
phase space for low energy scattering processes. The 
overlap between the quasi-particle state and the bare 
fermion state defines the wavefunction renormalization 
factor Z , which is suppressed from the non-interacting 
value of 1 but remains finite. Z is also the quasiparti¬ 
cle weight determining the discontinuity of the fermion 
occupation number at the Fermi surface. 

The interactions among quasi-particles are captured 
by the phenomenological Landau interaction function, 
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which describes the forward scattering processes of quasi¬ 
particles. The Landau interaction function can be de¬ 
composed into a set of Landau parameters F/ in which l 
denotes the partial wave channels. The physical observ¬ 
ables, such as compressibility, specific heat, and magnetic 
susceptibility, compared with their values in free Fermi 
gases, are renormalized by these Landau parameters. 

The Fermi surface can be made analogues to an elastic 
membrane. The energy cost to deform the Fermi surface 
can be viewed as the surface tension, which consists of 
two contributions from the kinetic energy and the inter¬ 
action energy. The kinetic energy cost is always positive, 
while the interaction part can be either positive or neg¬ 
ative. If the Landau parameter F[ is negative and large, 
i.e., Fi < —(21 + 1), then the surface tension vanishes 
in this channel, and then spontaneous distortion will de¬ 
velop on the Fermi surface^. This class of Fermi surface 
instability is denoted as Pomeranchuk instability in the 
literature. The simplest Pomeranchuk instability is fer¬ 
romagnetism which is an instability in the s-wave spin 
channel. 

The Landau interaction function also gives rise to col¬ 
lective excitations which are absent in free Fermi gases, 
such as the zero sound mode. The zero sound is a gen¬ 
eralization of the sound waves in fluids and solids. In 
fluids, the sound wave describes the propagation of the 
density vibration p(r,t), which is a scalar wave; in solids 
the sound wave is the vibration of the displacements of 
atoms from their equilibrium positions u(r,t), which is 
a vector wave. Compared to ordinary fluids which can 
only support density fluctuations, Fermi liquid possesses 
a microscopic structure of the Fermi surface which can 
be viewed as an elastic membrane, whose degree of free¬ 
dom is infinite described by the spherical tensor vari¬ 
ables Sni m . Consider a macroscopically small and micro¬ 
scopically large volume around r, around which a local 
Fermi surface can be defined. The local Fermi surface de¬ 
formation can vibrate and propagate, and thus generat¬ 
ing sound waves Sni m (r, t), which is the physical picture 
of the Fermi liquid collective excitations. The restoring 
force for the zero sound arises from Landau interactions 
rather than hydrodynamic collisions for the sound modes 
in ordinary fluids. 

B. Single-particle properties 

Let us neglect the influence of the confining trap, and 
also assume that dipoles polarize along the z-axis. The 
second quantized Hamiltonian of a single component elec¬ 
tric dipolar fermion system reads 

H d = Xl[e(fc)-Ii]c t (fc)c(fc) + ^- Y Vd{q) 

k k,k' ,q 

x ip\k + q)ip^ (k')tp(k' + q)ip(k). (12) 

In Sect. EH we define the dimensionless parameter 
as A = d 2 mkf / (3tt 2 H 2 ). It describes the interaction 





FIG. 2: The deformed Fermi surface of the 3D dipolar system 
at A = by the perturbative (solid) and variational (dashed 
red) approaches. The external electric field lies along the z- 
axis. From Ref. Q- 

strength, which equals the ratio between the average in¬ 
teraction energy and the Fermi energy up to a factor at 
the order of one. 

The Fermi surface structure of an electric dipolar 
fermion system is uniform but anisotropic. Intuitively, 
the inter-particle distance along the z-axis is shorter than 
that along x and y- axes because the dipolar interaction is 
attractive (repulsive) along the z-(xy) direction, respec¬ 
tively. Consequently, the Fermi surface will be approxi¬ 
mately a prolate ellipsoid, elongated along the z-axis and 
compressed in the equatorial cry-plane, which has been 
investigated in Refs. [3T| - f34l |. 

The above picture can be confirmed from the explicit 
calculation of the fermion self-energy at the Hartree-Fock 
level. The Hartree term vanishes because it involves the 
spatial average of the dipolar interaction. The anisotropy 
of the Fermi surface can be determined from the Fock 
term, while the latter also depends on the actual shape 
of the Fermi surface, thus they are coupled together and 
should be solved self-consistently. Nevertheless, at the 
leading order, we approximate the Fermi surface as a 
sphere with the radius in the free space as kf 0 , and then 
Yi HF (k) can be calculated analytically!! as 

£ HF (k) = -2\E kf0 P 2 (cosd k ) (13) 

k f 0 

where E kfo = and I 3D {x) = j|{3ir 2 + 8 - ^ + 

3(1 ~J ^ In | 11. In the two-component dipolar Fermi 
gases, the Hartree term still vanishes and the Fock term 
only exists for the intra-component interaction, thus the 
Hartree-Fock self-energy remains the same. 

The anisotropic Fermi surface distortion is determined 
by solving the equation of chemical potential y as 

eo {kf) + T, HF (k f ) = y(n,\), (14) 

where n is the particle density. The Fermi wave vector 
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kf depends on the polar angle as 

« 15 > 

in which the anisotropic distortion is at the linear order 
of A, and the A 2 term appears to conserve the particle 
numbers. Although Eq. [T5]is only valid at A <C 1, it pro¬ 
vides qualitative features. The Fermi surface anisotropy 
was also calculated by using the numerical variational 
method in Ref. [31]. The comparison between the ana¬ 
lytic perturbative result and the variational one is plotted 
in Fig. [2] for A = The Fermi surface based on the 
first order perturbation result Eq. [15] is less prolate than 
that based on the variational result. 

The anisotropy also exhibits in other single particle 
quantities. For example, the angular distribution of the 
density of states on the Fermi surface is calculated as 


Hartree-Fock level e k = e° + T,HF(k). The Landau in¬ 
teraction function is expressed at the Hartree-Fock level 
as 

f(k,k'-,q) = V(q)-V(k-k'), (19) 

where the first and second terms are the Hartree and 
Fock contributions, respectively. Due to the explicit 
anisotropy, /(fc, k'\ q) depends on directions of both k 
and k !, not just the relative angle between k and k' as in 
the isotropic Fermi liquids. 

The Landau interaction matrix elements for the dipo¬ 
lar system have been calculated in Ref. [33] by Fregoso 
et. al. According to the Wigner-Eckart theorem, the 
d r 2_ 3z 2 anisotropy of the dipolar interaction renders the 
following spherical harmonics decomposition as 


N(n k )^L = ^l^{i + ^\p 2 ( COS 0 k )}dn k , (16) 
47T n{ zirp 3 

where, N(£l k ) is the differential density of states along 
the direction of fl k . At the linear order of A, N(Q k ) 
develops the same anisotropy of P 2 { cos0). Thus the to¬ 
tal DOS at the Fermi energy does not change compared 
with that of the free Fermi gas. Nevertheless, it may be 
changed due to high order corrections. 


C. The Landau interaction matrix 

The anisotropic Fermi liquid theory has been con¬ 
structed in Ref. mum for the dipolar fermion sys¬ 
tems. The anisotropy of the interaction leads to the mix¬ 
ing among different partial-wave channels, thus we need 
to generalize the concept of Landau parameters into the 
Landau matrices. 

The variation of the Fermi distribution function at mo¬ 
mentum k is defined as 

dnj: = ng - n 0 ,k, (17) 

where no(k) = 1 — 9(k — k is the Fermi distribution 
function in the absence of interaction. The ground state 
energy variation of the single component dipolar Fermi 
gas is 


f(k, k'\ q) = ^2 


47T f w . r 


LI' 


\/(2i+l)(2i' + l) 


Ylm )Yl’m ) j 


( 20 ) 


where fw- m remains diagonal for the index m but couples 
partial wave channels with V = Z, l ± 2. The even and odd 
partial wave channels decouple because of the even parity 
of the dipolar interaction. The q dependence only appear 
in the channel of l = l 1 = m = 0, in which foo-,o(o) 

\", Other matrix elements at the Hartree-Fock level 
are tri-diagonal asiSS 


fifim = + (21) 

where 

(i) 47 t(Z 2 + l — 3to 2 )(2Z + 1) 


lm Z(Z + 1)(2Z + 3)(2Z — 1) ’ 

27 T i/[(Hri])T5r^2]|^TjI^jTTr^2j 


( 2 ) 

a L = — 


(Z + l)(Z + 2)(2Z + 3) 


• ( 22 ) 


For each l ^ 0, s satisfy the sum rule that 


SE : 


: 2 _ e k dn k 
k 


~rz~r mil ? 

2 Vo 9—>-0 

k,k' 


Y,fw=l;m = 0, (23) 


where, k, k' are momenta close to the Fermi surface; 
f(k, k!\ q) is the interaction function describing the for¬ 
ward scattering; q is the small momentum transfer for 
the forward scattering process, which is explicitly kept 
because of the non-analyticity of the Fourier component 
of Vd(q) as q —► 0; ng - = (cl + - c g)> which is reduced to 

the fermion occupation number as q —> 0; e k is the renor¬ 
malized anisotropic single particle spectra, and at the 


which reflects the fact the angular average of the dipolar 
interaction vanishes. 

To make the Landau matrix dimensionless, we multiply 
the single component density of states (DOS): = 

D fu?m, wh ere N$ D = (mk 3 °)/{2tnr 2 ) is the DOS 
of free Fermi gas and m* is the effective mass. At the 
linear order of A, m* = m. For concreteness, some low 
order Landau matrix elements are expressed at the linear 
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order of A as 


*oo;o(^V) — 
* 02;0 = 

* 22;0 = 
*22;±2 = 
*11;0 = 
*13;0 
*33;0 
*33;±2 


47tAP2(C0S 6q), 
—7tA; 

107r 


57T , 


J A, *22;±1 — — A, 

107T 

= -~ X; 

- 18 "x F - 9n \- 

- —A, ^ll;±l -- —A, 

= A, r 13i±I = 


5 
1477 


77T . 


— ■“r^A, *33;±i — — A; 


15 


= 0, 


10 

F , _ ?7r \ 

-*33;±3 —- —A. 

6 


(24) 


D. Thermodynamic quantities 

The thermodynamic properties, including the 
anisotropic effective mass and thermodynamic suscep¬ 
tibilities, are renormalized by the Landau interaction 
matrices. For simplicity, only the single-component 
dipolar systems are considered here. 


1. Anisotropic effective mass 


Landau parameters defined as follows 

*ii,n (0*) = 


*13,|| (^fc) = -7- E 


47T nt- *13;m 


=0,±1 V ^" 


Y* m (e k ,0)Y lm (0 k ,0). 


Thus to the linear order of A, the anisotropic radial ef¬ 
fective mass is 


1 1 — 7rAP2(cos0fe) 

m 3D,||( 0 *) 171 


(27) 


2. Thermodynamic susceptibilities 

Viewing the Fermi surface as an elastic membrane, we 
define the angular variation of the fermion distribution 
as 

/ 7.2,72, 

(2^3*"* (28) 

which can be further expanded in terms of the spherical 
harmonics as 

<5n(%) = E YimjQf^Snim. (29) 

lm 

For a Fermi surface distortion characterized by a set 
of Snim, the energy variation is calculated as 

= Me) H Sn *m Sn l'm K ll D ,m ~ hi m Sn lm j, (30) 

W;m 


It is well-known that the Galilean invariance leads to 
the relation between the bare mass of fermions and the 
effective mass of quasiparticles a a ZZ JS 

k_ = de{k)__ r d 3 k ' r ? 0n(e(fc')) 
m dk J (27r) 3 ^ } d P ' 

For an isotropic system, the effective mass m* is defined 
as de(k)/dk = k/m* for k on the Fermi surface. The 
renormalization of m*, or, the renormalization of the 
density of states, is = 1 + jjFf , which affects the 
specific heat as Cfl/Cfg = rn* /m, with Cfl and Cfg 
specific heats for the Fermi liquid and ideal Fermi gas, 
respectively. 

The dipolar Fermi gas is Galilean invariant so that Eq. 
E5l is still valid. However, due to the anisotropy, a self- 
consistent solution has to be done numerically. To the 
linear order of A, we approximate e(k') in the right-hand- 
side of Eq. [5S]with the free fermion energy. Defining the 

radial effective mass as rn* 3D ^(0 k ) = k (g k ) [k ■ we 

arrive at 


1 

m 


3D, 


(Ok) 


*11,11 (0fc) + *13,11 (Ok) , (26) 


where *n,|| (^fc) and *L 3 ,||( 0 fc) are the angular dependent 


where \o = —No is the density of states at the Fermi 
energy; fh* is the average value of the effective mass on 
the Fermi surface which equals m at the linear order of A; 
hi m = Mm + ^ 20 - Mm i s the external field in the partial 
wave channel of lm, and h\ 0 = |^/^A*j?^ is the explicit 
symmetry breaking arising from the dipolar interaction. 
The matrix kernel Km in Eq. [30] contains two parts as 


K u 


= Mi 


p3D 

r W -rr 


W: 


\J(21 + 1) (2Z' + 1) ’ 


(31) 


in which, Mw ]m is the kinetic energy contribution. At 
the linear order of A, Mw- m is calculated as 


M li?m = 2 + m\m S l,l’-2 


,( 1 ) 


( 2 ); 


S 2 )> 


where 


(i) _ I, ^ + 1 ) q (1) ,\ 

"Hm — i + a (m A\ a l.m.K 


4(2/+ 1) 


,(2) _ 3(Z + !)(/ + 2) 


(32) 


lm \\J(21 + 1)(2Z + 5) lr ' 

The expectation value of Sni m in the field of hi m can 
be straightforwardly calculated as 

Snim = Xo (Kf l ?. rn )- l h U m- (33) 


Thus Xo (KffiJ- 1 is the renormalized susceptibility ma¬ 
trix for the 3D dipolar Fermi system^!. 





















3. Pomeranchuk stabilities 


The inverse of an eigenvalue of the matrix Kw^ m can 
be considered as a thermodynamical susceptibility in the 
corresponding eigen-channel. If all the eigenvalues of 
Ku'-m are positive, i.e., this matrix is positive-definite, 
then the system is thermodynamically stable. If any of 
them becomes negative, the corresponding susceptibility 
diverges, which signals the onset of the Fermi surface in¬ 
stability of the Pomeranchuk typ o 34,43 . For the isotropic 
system, Ku>- m is already diagonal and to- independent, 
and the criterion for the Pomeranchuk instability is the 
well-known one 


Ft <-(21 + 1). (34) 

For the anisotropic dipolar system, the Ku'^ m matrix 
needs to be diagonalized to analyze its thermodynamic 
instabilities. The two strongest instabilities lie in the 
sectors of of to = 0 and to = ±2, and Z’s are even. For 
the case of to = 0, F 0 o- 0 (q) is singular as q —> 0, ex¬ 
plicitly depending on the direction of q. The most nega¬ 
tive eigenvalue occurs when q lies in the equatorial plane. 
The corresponding eigenvector mainly lies in the s-wave 
channel with 1 = 0, and the eigenvalue /i s reaches zero 
at A s = 0.135. This instability corresponds to the Fermi 
surface collapsing with a density modulation wave vector 
in the equatorial plane. This result nicely agrees with the 
numerical calculation in Ref. [3f|, in which the onset of 
an unstable collective mode implies the Fermi surface col¬ 
lapsing starting from A s ~ 0.14. It should be noted that, 
actually, this instability is mostly driven by the Hartree 
term interaction Vd(q), which cannot be simply dropped 
off by setting q = 0. Otherwise, the stability of the dipo¬ 
lar Fermi gas would be significantly overestimated. 

As for the sector of to = ±2, the eigenvectors of the 
minimal eigenvalues mainly lie in the d x 2 _ y 2 ± 2 i xy -w&ve 
channels. The lowest eigenvalues Pd± 2 touch zero at 
Xd ± 2 = 0.35. This instability corresponds to the biax¬ 
ial nematic instability of the Fermi surfaces studied in 
Ref. |43J. With the purely dipolar interaction, the s- 
wave channel instability occurs before the d- wave chan¬ 
nel one because Xd ± 2 > X s . Nevertheless, the s-wave 
channel instability can be cured by introducing a posi¬ 
tive non-dipolar short-range s-wave scattering potential 
Voo;Oj which adds to the Landau parameter of Fq^q with¬ 
out affecting other channels. 


E. The collective zero sound mode 


In this part, we review the calculation of the zero- 
sound-like collective mode in the dipolar Fermi gase s 34 ’ 36 . 
The anisotropic dipolar interaction brings a new feature: 
the zero sound excitation can only propagate within a 
certain range of directions beyond which the sound mode 
is damped. 


1. Generalized dynamical response functions 


The Boltzmann equation for the collective excitation 
of the single-component Fermi liquid isZS 


O 7* . 

—n(r, k,t) + ^ 


£ 


de(r, k , t) dn(r, k, t ) 
dki dri 

de(r , k, t ) dn(r , k , t) 
dri dh 


0, (35) 


where n[f, k, t ) and e(r, k, t) are the density and energy 
distributions in the phase space. 

In order to linearize the Boltzmann equation, the small 
variations of n(r, k, t ) and e(r, k, t) are defined as 

n(r, k, t) = n 0 ' d (k) + 6v(r,k,t), 

_ r (]3 ;,/ A A 

e( r,k,t ) = e H F(k) + J ^—^f(k,k')5v(r,k',t), (36) 

where 8v is defined with respect to the deformed equi¬ 
librium Fermi surface, and chf is the Hartree-Fock 
single particle spectrum. Substituting Su(r,k,t) = 
the linearized Boltzmann equation is 

arrived at 


5v(VL k ) 


N(Q k )Vk ' * 

uj-v,:-q 


dflk 1 

4-7T 


f m , 8 V(n kl ) = o, 


(37) 


where v? is the Fermi velocity; N(Q k ) is the differential 
density of states defined in Eq. [TGl 8v(Q, k ) is defined as 


5v(Sl k ) = J 8v(k ). 


(38) 


The spherical harmonics decomposition can be performed 
as 

^ T ^ ^ Xll"',mm'(tO, q)Fiffif- rn f'^8l / if rn / = 0, (39) 

I'm' l" 


where 8vi m is the component in terms of the spherical 
harmonics, and, 

XW,mm'(u,q) = ~ V{2l + m V+1) J dV p n , d P 

v 3D .q 

x Y? m (il p ) J W(IV). (40) 

to — v° u .q 

Due to the dipolar anisotropy, for a general propagation 
direction of q, , 5vi m in different channels are coupled. 


2. The s-wave channel approximation 

We first truncate Eq. SHI by only keeping the s-wave 
channel of l = 0. Even at this level, the anisotropy of 
the zero sound mode has already appeared. Taking into 
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8q/n 

FIG. 3: Dispersions of the zero sound s(9 q ) = Lo(9 q )/vf®q 
for the pure dipolar interaction at a) A = l/(37r 2 ) and b) 
A = 1/7T 2 . When the sound speed hits the particle hole con¬ 
tinuum, the sound becomes damped. These results are in 
good agreement with the numerical study in Ref. [36j. (From 
Ref. [Sf].) 


account the Hartree-Fock single particle spectra and the 
anisotropic Fermi surface, Xoo;Oo(w,g) is given by 


Xoo ; oo(w,<z) = 1- 


dflk N(Slk) s 

47 t No s - +ir)' 

(41) 


3. Correction from the coupling to the p-wave longitudinal 
channel 


Even in the isotropic Fermi liquid state, because the 
propagation direction q of the zero sound already breaks 
the 3D rotational symmetry to the uni-axial one, actually 
the zero sound mode mixes all the longitudinal channels 
of Snio- If the Landau parameter F\ is not small com¬ 
pared to Fo, the mixing between the s and p -wave longi¬ 
tudinal channels significantly modifies the sound velocity. 
In the isotropic Fermi liquid state, the modified sound 
velocity is determined by the following equation as^ 


-1 


Fo + 




1 + s 


1 — s ' 


(43) 


For example, in the 3 He system at 0.28 atm, if only con¬ 
sidering the s-wave channel, the sound velocity is cal¬ 
culated as s = Vg/vf = 2.0 based on Fff = 10.8. After 
including the coupling of the p -wave longitudinal channel 
in which F* = 6.3, the revised value of s increases to 3.6 
in agreement with the experimental measurements 79 . 

The case of the dipolar Fermi gas is more complicated. 
If the propagating direction q is not along the 2 -axis, 
no rotational symmetry is left, and thus, in principle, 
the longitudinal and transverse p -wave components are 
mixed. Here, the spherical harmonic functions Y) = i iTO= o 
( longitudinal) and Yi-i m -±i (transverse) are defined ac¬ 
cording to the principle axis along q instead of the 2 -axis. 
Nevertheless, usually the transverse p-wave channel mode 
is overdamped unless the p -wave channel Landau param¬ 
eter is positive and large, thus their effect to the zero 
sound mode is small and will be neglected. 

By keeping the mixing between the s-wave and the 
longitudinal p-wave modes, Eq. 0U1 is reduced to a 2 x 2 
matrix equation, and the collective mode can be solved 
based on 


Det[l + JV(g,w)] =0, 
where the matrix kernel of N (g, u>) reads 


(44) 


where = 9 v k^HF(k) , anc j propagation di- 

rection q is chosen in the ai 2 -plane with the polar angle 
8 q . Then the zero sound mode dispersion is determined 

by 


1 + -Fd0;0(Hg)X'00;00(w, (?) — 0. (42) 


The quantity s{9 q ) = ui(9 q )/(vf 0 q) is defined to repre¬ 
sent the angular dependent zero sound dispersion, which 
is solved numerically and plotted in Fig. [3] along with 
the edge of particle-hole continuum. The zero sound 
propagation angle is restricted and its dispersion s(9 q ) is 
anisotropic. For large angles of 9 q , the sound excitation 
enters the particle-hole continuum and is thus damped. 


N(q,u) 


Xoo ; oo (q, w)F 0 o ; o (q) 
XiO;00 (q, w)Foo;o ( q ) 


XlO;Oo(9,w) 

Xll;Oo((f,w) 


Fi 10 (g) 

_ 3 
Fii-o(q) 


(45) 


Firoiq) is the longitudinal p-wave Landau parameter de¬ 
fined as 


Fii;o{q) = cos 2 6> g Fii ;m= o +sin 2 0qFii ;m==F i, (46) 
and the response functions are 

dfl k N(fl k ) (q ■ k)f(Qk, fig) 


XiO;Oo(<?,w) = -VS 

Xll;0Cl((f, w) = -3 / 


' 47r No s-/(fifc,fi,) 
rffifc A^fifc) (q ■ k) 2 f (fife, fi q ) 
4tt N 0 s-f{ fifc,fi g ) 
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where q lies in the xz-plane with the polar angle 9 q , and 
q ■ k = sin 6 q sin 9 k cos <j>k + cos 9 q cos 9k ■ 

The numeric solution taking into account the 
anisotropic Fermi velocity and Fermi surface is per¬ 
formed, and the zero sound velocity as a function of 9 q is 
plotted in Fig. [3] The longitudinal p -wave mode modifies 
the sound velocity dispersion significantly. These results 
are in a good agreement with a fully numerical calcula¬ 
tion based on the same Boltzmann transport theory— 
This indicates that the zero sound mode is well captured 
by the coupling between the s-wave and longitudinal p- 
wave channels. 


IV. THE SO COUPLED FERMI LIQUID 
THEORY OF THE MAGNETIC DIPOLAR 
FERMIONS 

The magnetic dipolar interaction brings a new ingre¬ 
dient to the Fermi liquid theory, i.e., the SO coupled 
natur e 43 ' 58- — . Certainly, for the experimental system 
of the 161 Dy atoms whose hyperfine spin is as large as 
F = the theoretical analysis on magnetic dipo¬ 

lar interactions will be very challenging. Nevertheless, 
the spin-i case exhibits nearly all the qualitative fea¬ 
tures of the magnetic dipolar interaction, and thus will 
be used as a prototype model below. 

If the magnetic dipolar systems are partially polarized, 
this SO coupling already appears at the single-particle 
level exhibiting anisotropic Fermi surfaces, which is a re¬ 
sult from the ferro-nematic coupling as shown by Fregoso 
et. a/ i 43 i 59 . For the unpolarized case, the Fermi sur¬ 
faces remain spherical without splitting. Nevertheless, 
the effects of the SO coupling appear at the interaction 
level 61 ' 62 , including Fermi surface Pomeranchuk instabil¬ 
ities, and SO coupled collective modes. 

The second quantized Hamiltonian of the spin-A 
fermions with the magnetic dipolar interaction is ex¬ 
pressed as 

Hmd, = ^ ~ A t ] c a(fc)c Q .(fc) + r 'y ' F a p-p'ot'{q) 

k,cx. k,k' ,q 

X iplik + qjiplik’^p'ik' +q)-ip a '(k). (47) 

Similarly, the dimensionless interaction parameter can be 
defined accordingly as A m = p^gpmkf 

The SO coupled Landau interaction function is re¬ 
viewed and the SO coupled partial wave decomposition 
is performed in Sect. IIV Al the Pomeranchuk instability 
is reviewed in Sect. IIV Bl the zero sound-like excitation 
with the SO coupled feature is reviewed in Sect. IIV Cl 


A. The SO coupled Landau interaction 

In Ref. [61 j, the Landau interaction function of the 
magnetic dipolar fermions with a general hyperfme-spin 


F was given. For simplicity and concreteness, below we 
still use the spin-i case for illustration— Based on the 
Fourier transform of the magnetic dipolar interaction Eq. 
OH the Landau function of the spin- A magnetic dipolar 
system is expressed at the Hartree-Fock level as 


faa',pp'(k,k’-,q) — f aa ' t pp'(q) + faa',pp'(k, K ) 


7rg 2 F n 2 B 


{^MoLOt' ,PP' isi) (p)) } 


where k and k! are at the Fermi surface; the small mo¬ 
mentum transfer q in the Hartree term is explicitly kept 
due to the singularity at q —> 0 in Eq. [11] p is the 

unit vector defined as p = , ; the matrix kernel 

\k—k'\ 

M ac only depends on the direction of rh as 


M aa i t pp>(rh) = 3 {3 aa , ■ rh)(app> • rh) - ff aa > ■ app>, (48) 

for rh = p and q. In order to arrive at f£ a , pp,(k, k') in 
Eq. HU the following identity is used 

‘ p)(V3 oi' * V) @ap' ' & pa' 

= 3(dw • p){vpp’ • P) - • app>- (49) 


1. The SO partial-wave decomposition 


It is convenient to work in the SO coupled bases for 
the magnetic dipolar Fermi liquid theory. The varia¬ 
tion of the single particle density matrix in momentum 
space is defined as Sn aa >(k) = n aa f(k) — 8 aa 'no(k), where 
n aa ’(k) = (?/4,(fc)^ a '(fc)) and no(fc) refers to the ground 
state distribution of the free Fermi system. As for spin 
indices, Sn aa ' (k) can be expanded as 

Sn aa > ( k) = ^2 Sn Ss z (k)xss z , (50) 

Ss z 


where xss z ,aa’ are the bases for the particle-hole singlet 
(density) channel with S' = 0 and triplet (spin) channel 
with S = 1, respectively, defined as 


Xl0,aa / 


\00,aa / 


®z,aa'i Al±l ,aa' 
8 aa '. 


T 1 / 
V5 ( ' 


± ia 


y,czoi') 


(51) 


Similarly as before, we integrate Sn aa '(k) along the 
radial direction, and arrive the angular distribution 
5n aa '(k). In the SO decoupled bases, it is expanded as 


8riaa' (*)= £ 8riLmSs z 8i Lm(k')XSs z ,aa' ■ ( 52 ) 

LmSs z 


More conveniently, it can be reorganized in the SO cou¬ 
pled bases as 


Sn 

aot' (*) = £ 8 tijj z -ls yjj z -,Ls(k,aa'), (53) 

JJ Z ;LS 
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where 6njj z;LS = J2ms z ( LmSs z\ J 
yjJz;Ls(k, aa') is the SO coupled spherical harmonic 
functions 


y’j Jz;Ls(.ki Oi(X ) ^ ' (LlTl'Ss z \ J J zf^Lm (k^XSsz ,aa' • 

msz 

Using the SO coupled bases, the Landau matrix is di¬ 
agonal with respect to the total angular momentum J 


and its ^-component J z as 


N, 


-^-faa'-,pp'(,k,k') = ^2 yjJz-Li{k,aa')Fjj zL i.jj zL 'i 


J J Z LL' 


x yjjz-,L'i(ki PP')- 


(54) 


The matrix kernel Fjj z j^jj z jj\ reads 


F 


ttX 


JJzLl-JJzL'l — 


6j,i6 l ,oS L ',o(.26j Z! o-Sj z , ±1 )+ ^ (Lmls z \ JJ z )(L'm'ls' z \ (55) 


in which the first term is the Hartree contribution with q set as the 2 -axis; the second term is the Fock contribution 
with T[ rnls ^. L , m , ls , defined in the SO decoupled bases as^>^ 


rjiF 

Lmls Z ]L'm' ls' z 


6l- 


2 ,L' 


ttA/ 5ll> _ 5l+2,l' 

2~\L(L + 1) ~~ 3(L + l)(L + 2) 3(L - 1)L 

x Y Lm (n r )YZ, m ,(n r ). 


) J dn r (5»,y,-47ryi 8 ,(n r )y 1 * i (n r )) 


(56) 


The same value of J may arise from L with differ¬ 
ent parities. Below we use ( J ± J Z LS ) to represent dif¬ 
ferent angular momentum channels, where ± is the par¬ 
ity eigenvalue. The Hartree term of Eq. [55] only con¬ 
tributes to the (1 + J z 01) sector, which explicitly depends 
on J z because q (chosen as the 2 -axis) breaks the 3D 
rotation symmetry down to the uniaxial rotation sym¬ 
metry. For other sectors, Fj±j zL1;J ± j z l’i does not de¬ 
pend on J z as required by the Wigner-Eckart theorem. 
The matrix Fjj z li-jj z l'i is nearly diagonalized except 
for the case with L ^ L'. For J + = 0 + ,2 + ,.. and 
J~ = 1 — , 3 ,..., there is only one possibility that L = J , 
and thus Fj± j z li,j±j z l ’1 is already diagonalized. In 
comparison, for J + = l + ,3 + ,... and J~ = 0 - ,2 - ,... L 
can take two different values of L = J ± 1, and thus the 
matrix is reduced to 2 x 2 (For J~ = 0“ the only possi¬ 
bility is L = 1.). The concrete forms of Fj± nrw_± j z l'\ 
for a few low orders of are given in Ref. [61 II 621 . 

2. Thermodynamics susceptibilities 


Xo = No is the Fermi liquid density of states; hjj z LS is 
the external field. At the Hartree-Fock level, N 0 receives 
no renormalization from the magnetic dipolar interac¬ 
tion. The expectation value of Stijj z ls is calculated as 

Sn JJzLS = Xo '52( M )jJzLS;JJzL'S h JJ*L'S- ( 59 ) 
L' 

The thermodynamic stability condition is equivalent to 
that all the eigenvalues of the matrix Mjj z ls-jj z l'S are 
positive. 

Fregoso et a h 43 i 59 found that a uniform magnetic field 
h along the 2 -axis not only induces spin polarization, but 
also a spin-nematic order in the channel of (J + J Z LS) = 
(1 + 021). As a result, the external magnetic field induces 
an effective SO coupling 

Hhso = ^Tr\h^2'ipi(k){[(k 2 -3k 2 z )a z 
k 

- 3k z (k x a x + k y a y )\ }i/jp(k). (60) 


The variation of the ground state energy in the SO 
coupled bases is expressed as 


SE 

~V 


167r| 

E 


E Sn JJz;LS M JJzlS-,JJzL'sSn JJz 
JJ Z LL'S 

hjj z Ls5njj z -LS f, 


;L'S 

(57) 


JJzLS 

where the matrix kernel is 

Mjj z ls\JJzL's = 5lu + Fjj zL s-jj z l's\ (58) 


Apparently, Eq. [201 breaks TR symmetry, which is 
markedly different from the relativistic SO coupling in 
solids. 


B. SO coupled Pomeranchuk instabilities 

If one of the eigenvalue of the Landau matrix becomes 
negative, Pomeranchuk instability occurs in the corre¬ 
sponding channel. The SO coupled nature in the mag¬ 
netic dipolar system manifests as follows. 
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Let us first consider the channel of J = 1 + — . In the 
absence of external fields, a density wave with a long 
wave length q —> 0 can take the advantage of the Hartree 
channel interaction. The leading instability lies in the 
sector with J z = ±1. Without loss of generality, we use 
the J z = 1 sector as an example: 


( Aiioqiioi 

Elioi;1121 ' 

i 7rA m j 

^-4v / 2\ 

V Eli21;1101 

Eli21;1121 J 

1 12 ' 

1 J 


whose negative eigenvalue and the associated eigenvector 
are 


w\ +1 = -0.377rA m , ipl +1 = (0.97,-0.25) t (62) 

The instability occurs at u>( +1 < —1, or, equivalently, 
A m > AJ +1 = 0.86, and the eigenvector shows that it is 
nearly a ferromagnetic instability hybridized with a small 
component of the ferro-nematic channel. The spin polar¬ 
izations lie in the xy- plane and the spiral wave vector 
is along the z-axis. The spiral wave vector q should be 
at the order of the inverse of the system size in order 
to minimize the energy cost of spin twisting, which may 
further depend on concrete boundary conditions. 

An interesting point is that because of the coupling 
between the ferromagnetic and ferro-nematic orders, the 
Fermi surfaces are distorted^. The oscillations of the 
distorted Fermi surfaces couple to spin waves and bring 
Landau damping to spin waves. This may result in non- 
Fermi liquid behavior to fermion excitations. Similar ef¬ 
fects in the nematic Fermi liquid state have been exten¬ 
sively studied before in the literature^ - — . 

The next sub-leading instability is in the J = 1~ chan¬ 
nel with L = 1 and <5=1, which has been studied in Ref. 
pm. For J z = 0, the generated SO coupling at the sin¬ 
gle particle level exhibits the 3D Rashba-type as 

H so ,i- = \'^z\y^ J ‘>Pa(k)(k x ay - k y <j x ) a pipp(k), (63) 

k 

where \n z \ is the magnitude of the SO order parameter. 
The magnetic dipolar interaction already possesses the 
SO nature. The Pomeranchuk instability brings it to the 
single particle level by breaking rotational symmetry and 
parity. 

The instability of the J = 1“ sector is similar to but 
fundamentally different from the Pomeranchuk instabil¬ 
ity in the L = S = 1 channel studied in Refs. [Hill- 
Different from the magnetic dipolar systems, the starting 
Hamiltonian in the latter case has no SO coupling at all. 
After instability occurs, an effective SO coupling appears 
at the single particle level. In particular, a k ■ a type 
SO coupling can be generated with total angular mo¬ 
mentum J = 0. It is generated from the SO decoupled 
interactions through a phase transition, and thus was 
denoted as spontaneous generation of SO couplin g 84 ' 85 . 
They are a particle-hole channel analogy of the p- wave 
triplet Cooper pairings of the isotropic 3 He-I? phased 
and the anisotropic 3 He-A-phas o 87 ' 88 , respectively. 


C. The SO coupled collective modes — the 
topological zero sound 


We review the study of the collective modes in the 
SO coupled Fermi liquid theor y 61 ' 62 . The generalized 
Boltzmann equation including the spin degree of freedom 
can be written as2i 

—n(r,k,t ) - ~[e(r,k,t),n(r,k,t)] 

1 ^ r de(r , k, t) dn(r, k , t) t 

2 l dki ’ dVi J 

i 

1 r de(f, k , t) dn(r , k , t) 

2 l dri ’ dki 

i 



where n aa ' (r, k, t) and e aa < (r , k, t) are the density and en¬ 
ergy matrices for the coordinate (r, k ) in the phase space; 
[,] and {,} mean the commutator and anti-commutator, 
respectively. After linearizing Eq. [BJ] and expanding it 
in the plane-wave bases 


fin aa > (r, k) = 5n aa ' (k)e iq ' r 


(65) 


we arrive at 

- 1 cos 9k f , r N 0 

dn aa '(k) — - - „ y y j dllk'—Jaa',t 


2s — cos dk “ 
x 6npp>(k')= 0, 


'(MO 

( 66 ) 


where s is the dimensionless parameter w/ ( v/q ). Without 
loss of generality, the propagation direction of the wave 
vector q is defined as the z-direction. 

By using the SO coupled bases &njj z -LSi Eq. [111 is 
rewritten as 


Snjj z -LS + ^2 KjJzLS,J'J x L's{s)Fj'j z L'S-,J'J z L"S 

J';LL' 

X Snj' Jz L"s = o, (67) 

The matrix kernel Kjj z ls-,j'j z l'S reads 
Kjj z lS;J'j z l's{s) = (LmSs z \ J J z ) (L'mSs z \ J' J z ) 

ms z 

X D LL / ;m (s), (68) 

where £Ill’(s) is equivalent to the particle-hole channel 
Fermi bubble in the diagrammatic method defined as 

«LL';m(s) = - / dn k Y£ m (k)Y L , m (k) ^ \ ■ (69) 
J S - cos 6 k 

The largest positive Landau parameter lies in the 
(1 + 001) channel, which can support propagating modes. 
Since q breaks parity and the 3D rotation symmetries, the 
(1 + 001) channel couples to other channels with J z = 0. 
The Landau parameters of orbital partial wave channels 
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0 k x /k F 1 


FIG. 4: The spin configuration (Eq. 1721) of the zero sound 
mode over the Fermi surface shows the hedgehog type topol¬ 
ogy. Although the hedgehog configuration is distorted in the 
z-component, its topology does not change for any values of 
A m describing the interaction strength. From Ref. [89]. 


and n 0 - 0 n(q) = J2k + <t)(& ' ^typify, the former 
mode describes spin oscillation along the direction of q ,, 
and the latter exhibits a hedgehog configuration of spin 
distribution on the Fermi surface. The eigen-mode is a 
hybridization between them, which can be represented as 

( u 2 sin 6G cos ^ \ 

U 2 sin Ofe sin <f>% e i&?-«iv f t), ( 72 ) 

U 2 cos + Ui ) 

where (ui,u 2 ) T is the eigenvector for the collective mode. 
For all the values of A m , \u 2 \ > |iti| is satisfied. Thus, the 
spin configuration, as shown in Fig. [IJ is topologically 
non-trivial with the Pontryagin index ±1. The sign of 
the Pontryagin index periodically flips at the nodes of 
the sound wave as the time and spatial coordinates varies. 
This collective mode can be considered as the topological 
zero sound. 


with L > 2 are small, which will be neglected henceforth. 
There are three spin channel modes with L = S = 1 and 
J z = 0 denoted as (0 — Oil), (1 — Oil), and (2“011). Even 
in the presence of q ., the system still possesses the reflec¬ 
tion symmetry with respect to any plane including q. By 
writing down their bilinear fermion expressions, we can 
check that among the above modes, (1 + 001), (0 _ 011), 
(2 _ 011) are odd, and (1”011) is even under this reflec¬ 
tion operation. Furthermore, the (2 _ 011) channel can be 
neglected because the Landau parameter in this channel 
is about one order smaller than those in (1+001) and 

(o-on). 

Now, we only keep two coupled modes (1+001) and 
(0”011). Using the following relations 


Aiooi ; iooi(s) 

-Kl001;001l(s) 

Aooii.oon(s) 


^00;o(s), 

A0011;100l(s) = S^00;C)(s), 

I(lml - m|00)| 2 Un ;m (s) 

m 

^00;o(s)> (70) 


we reduce the coupled 2x2 matrix equation based on 
Eq. [67] into 

^oo ; o0) = 1 - ^ln|^-t-^| + j^s0(s < 1) 

2 1 — s 2 

Ft ± Jf$ + 4(s 2 - 1)F* 

- V- i )F> —■ < 71 > 


where 

F + = T'looiiiooi + ^0011;00111 = ^1001:1001-fboil;0011- 

For the two branches of Eq. [TTJ only the one with the 
minus sign possesses the solution with s > 1, as required 
by the condition of the undamped collective mode of the 
Fermi liquid. Since n 1+001 (g) = 


V. UNCONVENTIONAL TRIPLET COOPER 
PAIRING WITH MULTI-COMPONENT 
DIPOLAR FERMI GASES 


The p -wave (L = 1) triplet (5 = 1) Cooper pairing is a 
celebrated unconventional pairing superfluid state which 
has been a research focus of condensed matter physics 
for decades (see Ref. (5filf9fli | for reviews). The typi¬ 
cal system is the superfluid 3 He which exhibits both the 
isotropic B phased and the anisotropic A-phas o 87 ' 88 . So 
far, the most accepted pairing mechanism is spin fluctua¬ 
tions arising from the prominent ferromagnetic tendency 
because there exists a hard core part in the interaction 
between two 3 He atoms. 

As explained in Sect. [U the anisotropic dipolar in¬ 
teraction provides a novel and robust mechanism for 
the p -wave spin triplet Cooper pairing and the compe¬ 
tition between the singlet and triplet channel pairing 
instabilitie o 28 ’ 29 ! 44 - 54 ! 69 . Furthermore, the coexistence of 
the singlet and triplet pairing symmetries in the electric 
dipolar systems naturally leads to a novel TR symmetry 
breaking mechanism first pointed out in Ref. |50(. Re¬ 
cently, this mechanism is also studied in a wide context of 
superconducting system s — * 92 . For example, it has been 
later proposed in the iron-based superconductors^, and 
its topological electro-magnetic responses have been also 
studied^. In fact, the TR symmetry breaking pairing is 
also a research focus in condensed matter physics, such as 
in the study of high T c cuprates (e.g. YBa2Cu306+a;)^ 
and ruthenates (e.g. Sr 2 Ru 04 ) 2 i. 

In this section, we will first present a brief overview of 
the p -wave pairing in the superfluid 3 He system in Sect. 
IV Al The pairing symmetries in the single-component 
and two-component electric dipolar fermion systems are 
reviewed in Sect. IV Bl and Sect. IV Cl respectively. The 
effect of the TR reversal symmetry breaking is reviewed 
in Sect. IV PI 
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A. A quick overview of the p-wave triplet pairing 
in the superfluid 3 He system 

The general structure of the p -wave spin triplet pairing 
takes the form of a 2 x 2 symmetric matrix A a p(k) in 
momentum space, which can be conveniently represented 
by the d-vector as^ 


Yi(k-k’) V d (k+k’) 



A a p(k) = A(A-)d Al (fc)(icr M cr y ) Q , /3 , (73) 

where A(fc) is a complex number; d describes the spin 
degree of freedom, which is a complex unit 3-vector sat¬ 
isfying d* ■ d = 1; ia^cr v (pL = x,y,z) form the bases of 
the 2x2 symmetric matrices for the triplet pairing. The 
so-called unitary pairing means that A^(k)A(k) equals 
the identity matrix up to a constant, i.e., A (k) is pro¬ 
portional to a unitary matrix. In terms of the d-vector 
language, the unitary pairing corresponds to d* x d = 0, 
i.e., d is equivalent to a unit real vector up to an overall 
complex phase. 

The two most prominent superfluid phases of 3 He are 
B and A phases, both of which belong to the class of uni¬ 
tary pairing. In the B phased, a typical pairing matrix 
structure is 

A^(fc) = Ak ' (°iv v )a f 3, (74) 

where A is a complex constant. In this configuration, L 
and S are combined into J = 0, and thus, the B phase 
is rotationally invariant and also a fully gapped phase. 
In other words, the B phase spontaneously breaks the 
relative SO symmetry, exhibiting the SO coupled paired 
structure. It maintains TR and parity symmetries. Fur¬ 
thermore, it has also been recently found that the B 
phase is actually a 3D topological superfluid phase. In 
comparison, the 3 He-A phase is anisotropic, in which L 
and S are decoupled and J is not well-define d 87 ’ 88 . A 
typical pairing matrix for the A phase is 

A^(k) = A(k x + ik y ){d ■ aia v ) a p, (75) 

in which d is an arbitrary unit 3-vector. It breaks TR 
symmetry and exhibits nodal points on the north and 
south poles on the Fermi surface. 

We will see in this section and in Sect. IVIl that both the 
electric and magnetic dipolar fermion systems support 
novel p- wave triplet structures which are different from 
the 3 He-I? and A phases. 


B. The pz-wave pairing with the single-component 
dipolar fermions 


FIG. 5: Pairing interaction matrix elements for the single 
component dipolar interaction V t (fc; k') = Vd(k — k') — V d (k + 
k') due to the two-particle interference and the Fermi statis¬ 
tics. 

In real space, this p 2 -pairing symmetry is also clear 
since the electric dipolar interaction is the most attrac¬ 
tive if the displacement vector between two fermions is 
along the z-axis. Below, we present the partial wave anal¬ 
ysis in momentum space. The pairing interaction can be 
expressed at the mean-field level a o 28 ’ 29 

H pa i r = -^r^2v d (k-k')^(k)^(-k)ip(-k')i/j(k'), 

2V ° k ,w 

= A TtH (76) 

where 

V t (k;k') = V d (k-k')-V d (k + k') (77) 

satisfying V t (k',k') = V t (—k;k') = —V t (k,k') as shown 
in Fig. [5] The symbol V t is used because the pairing 
analysis here is the same as that for the triplet pairing of 
the two-component dipolar fermions in Sect. IV Cl 
The BCS mean-field gap equation is 

_ r /jSi.' _ i 

m = - j J^fV t {k-k')[K{k')-—\ A(k'), (78) 

where E(k) = \J (e* — p) 2 + A 2 (k) and K(k) = 

tanh.[^E(k)]/[2E(k)]. The integral of Eq. [75] has been 
regularized following the f-matrix method at the level of 
the Born approximation. This regularization is equiva¬ 
lent to truncate the energy away from the Fermi energy 
±d>, and ui is at the order of the Fermi energy^. 

Although the p z -pairing is very intuitive in the real 
space, it is not so obvious in momentum space. Before 
performing the partial-wave analysis, a qualitative mo¬ 
mentum space picture can reveal why the p 2 -pairing is 
natural as shown in Fig. [G] (a). Let us set k || z and 
k' —> k, then (k — k') L z and (k + k') || z, thus 


The pairing symmetry structure of the single¬ 
component dipolar system was studied in early works of 
Refs. [H|j29[. In this case, the gap function is simplified 
as a complex number, and thus, the d-vector notation is 
not needed. 


V t (k; k') = -4ird 2 < 0, V t {k\ —k') = And 2 > 0, (79) 

which favors the pairing in polar regions with an odd 
parity. On the other hand, if we set k in the equatorial 
plane, say, k || x. and also k! —> k , then k — k! lies in 
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FIG. 6 : Pairing interaction matrix elements for the dipolar 
interaction Vt(k;k'). a) For the case of k || z, and k' —> k, 
Vt(k\k') = Vt(k\~k') < 0. b ) For the case of fc || x , and 
k' —»• k, Vt(k- k') varies and its angular average is positive. 


the yz-plane. The value of V d {k — k') depends on the 
polar angle of the vector of k — k!. Its average is ^ d 2 
and that of V d (k + (?) = -^-d 2 , thus the average of 

Vt(k ; k') as k' — > k for k in the equatorial plane is positive. 
This means that the pairing amplitude is suppressed in 
the equatorial plane. Combining the pairing structures 
both in the polar region and the equatorial plane, it is 
clear that the pairing symmetry is mostly of the p z - wave, 
which is also consistent with the real space picture of 
pairing. 

Around T c , the gap equation can be linearized. The 
standard pairing eigenvalue analysis is performed by 
defining the eigenvectors (\>\ s (k) satisfying the eigen- 
equation 

N 0 J dn k ,V t (k; k'Wk') = w\ 4 iS (k), ( 80 ) 


where w\ are eigenvalues; _/V 0 is the density of states on 
the Fermi surface; i is the index of eigenvectors. The 
actual pairing occurs in the channel with the lowest neg¬ 
ative eigenvalue. The spherical harmonics decomposition 
of Vt jS (fc; k') reads 

^v t (k;k')= J 2 ( 81 ) 

1,1'\m 


where l and l' only take odd values. Compared with 
Fa nTJl and Eq. [20J we arrive at 


Vll'\m 


___2iVo fu/ym__ 

\J (21 + 1)(2 V + 1) 


(82) 


The diagonalization shows that the most negative 
eigenvalue lies in the sector with m = 0 , whose eigen¬ 
value and eigenvector are 


wi =° =0 = -3.82A, <(> 2 (flfc) » 0.99Yi 0 - 0.12y 3O . (83) 

The above pairing eigenvector shows that the pairing 
symmetry is mainly of the p z - type, in agreement with the 
intuitive real space analysis. The pairing eigen-values in 
the sector of to = 0 can even been solved analytically as 
shown in Ref. [29| as w l t rn=0 = — 12 A/[ 7 r( 2 i + l) 2 ]. The 


gap function of Eq. [83] vanishes for quasi-particle mo¬ 
menta lying in the equatorial plane, and thus exhibits a 
nodal plane in the Bogoliubov excitation spectra. 

As shown in Ref. [29[, the standard mean-field value 
of T c is related to the eigenvalue solved above as 


2e 7 w ^ 
T c = -e“* , 


(84) 


where 7 = 0.5772 is the Euler constant. In Ref. [29j . 
a further analysis based on Gor’kov, Melik-Barkhudarov 
(GM) approach*^ is performed, which takes into account 
the media polarization effect due to the virtual process 
of creating particle-hole excitations. This GM approach 
shows that the energy cutoff effect can be approximated 
by ui ~ 0.42e/. 


C. Competition between the triplet and singlet 
pairings in the two-component dipolar Fermi gases 


The new ingredient of the two-component dipolar 
fermions is that both the spin singlet and triplet pair¬ 
ings are allowed^- 3 -. We define the pairing operators in 
these two sectors as 


P s (k) 

P?$) 


1 

—ptr 

V 2 

1 

—ptr 

V 2 


P{k)(-ia v ) 


P{k)(—ia v cr^) 


(85) 


with p = x,y,z, respectively, where P a 0 (k) = 
it> a {k)tpp(-k). Then the pairing Hamiltonian is ex¬ 
pressed as 

H pair = E { [ E P^(k)P?(k') } 

+ V s (k-,k')Pt(k)P s (k') }, ( 86 ) 


where V t}S {k; k') = Vd{k—k')^fVd{k+k') are pairing inter¬ 
actions in the singlet and triplet channels, respectively. 
The Bogoliubov quasiparticle spectra become Ei(k) = 

\J ( e k — pi ) 2 + Af (fc), and A 2 ^{k) are the eigenvalues of 

the positive-definite Hermitian matrix A^(k)A(k). 

The gap equation takes the matrix form as 
_ r fj'iy _ 

A a/} (k) = -J j^3 V d (k-k , ){\'i/j a (k)'i/jp(-k)\), (87) 


where (||) means the thermal ensemble average. A a p(k) 
can be decomposed into the singlet and triplet channel 
pairings as 

A a p(k) = A s (k)i(J v a p + A t ^(k)(icr^<j v ) a p, ( 88 ) 

in which A s and A t, P satisfy 


A t(s),/j.(k) 


- 5 E / 

x~]A t(s),/j,(k ), 


( 89 ) 
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where Ki ( k ) = ta,nh[^Ei(k)\/[2Ei(k)\. After linearizing 
the gap equation around T c , we perform the eigenvalue 
analysis for the pairing problem. The eigen-equation of 
the triplet sector is the same as Eq. [5Ql while that of 
the singlet sector can be obtained by replacing V t {k: k') 
with V s (k: k'). The spherical harmonics decomposition 
of V s (k: k') can be done in a similar way to Eq. [HU and 
the resultant Vw-m takes the same form as Eq. [HU but 
with 1,1' only taking even integer values. However, the 
Hartree interaction does not exist in the pairing channel, 
and thus, for l = V = m = 0 , Voo ; o = 0 for the case of 
the purely dipolar interaction. 

The analysis of pairing eigenvalues in the triplet sector 
is the same as that in the single component case. Thus, 
the leading pairing symmetry still lies in the Pz-channel. 
The pairing eigenvalue and eigenvector are the same as 
Eq. 1831 and T c is still approximately determined by Eq. 
l84l We can express the triplet pairing in terms of the 
d-vector as 

A *,„(£) = A (90) 

where 7 is the U(l) phase and A t is the pairing am¬ 
plitude. This spin-triplet pairing breaks the U(l) gauge 
symmetry in the charge channel and the SU(2) symmetry 
in the spin channel, and thus, there are two different low 
energy excitations: phonon and spin-wave modes. This 
pairing is still invariant under a combined ^-symmetry 
of —» — d M and 7 —>■ 7 + and thus, it supports 

two different classes of vortices: the usual integer vor¬ 
tex of superfluidity, and the half-integer quantum vortex 
of superfluidity combined with a 7r-disclination of the d- 
vector. 

In the singlet channel, the lowest eigenvalue and the 
corresponding eigenvector are 

w° s = -1.93A, r +d (Vk) « O. 6 F 00 - O. 8 Y 2 O) (91) 

respectively. The eigenvector mixes the s and d fc 2 _ 3 fc 2 - 
channels. All other negative eigenvalues are small and 
negligible. Although for the purely dipolar interac¬ 
tion, the singlet channel pairing instability is significantly 
weaker than that of the triplet channel. Nevertheless, the 
matrix element Voo ; o may receive additional contributions 
from the short-range s-wave scattering interaction, which 
in principle is tunable through Feshbach resonances. At 
Voo ; o/A ~ —3.15, the singlet and triplet channel instabil¬ 
ities become degenerate. 

Considering the competition between the spin triplet 
and singlet channel pairings, we generally expect two 
pairing superfluid transitions. The first transition is trig¬ 
gered by the stronger pairing channel, say, the triplet p z - 
channel, which determines T c \. In the case that the sin¬ 
glet channel pairing is weaker but nearly degenerate with 
the triplet one, a second transition may occur at T c 2 as 
further lowering temperature. The mixing between the 
single and triplet pairings breaks parity, and thus, the 
second transition is also a genuine phase transition. The 


coupling between these two pairing channels can be cap¬ 
tured by the following Ginzburg-Landau (GL) free energy 
as £2 

A F = 7l (A? • A t )|A s | 2 + 7 2 {A: • A*A,A S + ccj 
+ 73 |A( x A t | 2 , (92) 

where A t is a compact notation for Aj iM (/j, = x,y,z); 
other non-gradient terms in the GL free energy only de¬ 
pend on the magnitudes of order parameters. 73 should 
be positive, such that A t can be described by a real d- 
vector multiplied by a U(l) phase. The sign of 72 deter¬ 
mines the relative phase between A s and A t : if 72 < 0, 
the phase difference between A s and A t is 0 or it; if 
72 > 0, the phase difference is ±f. The latter case gives 
rise to a novel pairing with TR symmetry breaking as 
explained in Sect. IV Dl 

In SO coupled systems, the coupling between the sin¬ 
glet and triplet pairings through a spatial gradient is 
considered in Ref. j98|, which leads to a spatially non- 
uniform state. Nevertheless, because of the spin con¬ 
servation, such a term is not allowed in electric dipolar 
systems. 

D. TR symmetry breaking mixing between the 
singlet and triplet pairings 

When the singlet and triplet channel pairings coexist, 
a natural question is: what is the relative phase between 
these two pairing channels? It was found in Ref. [50j 
that, at the mean-field level, the phase difference of 
between A s and A t,Al is favored, so that TR symmetry is 
spontaneously broken. This is a general mechanism lead¬ 
ing to the TR symmetry breaking pairing, which applies 
for the mixing problem between singlet and triplet pair¬ 
ings in the absence of SO coupling in the weak coupling 
limit. It has also been found later in Ref. [52ll53j | that 
the TR symmetry breaking also exists in the resonance 
interaction regime of the dipolar Fermi gases. 

The reason for the above TR symmetry breaking ef¬ 
fect is that the weak coupling theory favors the unitary 
pairing, i.e., A^(/c)A (k) is an identity matrix up to a 
factor. A simple calculation shows that A^(k)A(k) = 
|A S | 2 + |A ^| 2 + Re(A s ’*A t ’ M )(a M ) r . The unitary pair¬ 
ing is achieved if and only if a phase difference e 1 ^ = ±i 
exists between A s and A t,Al . More precisely, it can be 
proved following the method presented in Ref. [99|. The 
quasi-particle spectra read Ei = y/£ 2 + A 2 (i = 1, 2) with 

Ai j2 = |A s <(> s+d (H fe )| 2 + |A t <f(H fe )| 2 

± 2Re(A*A t )</> s+d (Hfc)^(Hfc). (93) 

The quasiparticle contributi on to the free energy 
is f(x ) = — j ln [2 cosh | + x\ , which satisfies 

> m. Thus, 

/(A 2 ) + /(A 2 )>2/(^±^t) (94) 
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and the minimum is reached at A^ = Aj. This is precisely 
the condition of the unitary pairing. 

Here is another intuitive view to see why the phase 
difference e 1 ^ = ±* is favored. Without loss of generality, 
let us set d || z, and assume a relative phase difference <fi. 
The ratio between the pairing amplitudes is 

(V> t (fc)^(-fc)) = A s + e^A t 
(ipl(k)^(-k)) t 

Only when e u ^ the magnitude of this ratio is 1, such 

that the pairing strengths are the same. Otherwise, say, 
if e l(t> = ±1, the pairing amplitudes of and 

(il>l(k)ijj^(—k)) are not equal. Again, the pairing struc¬ 
ture, whose gap magnitudes distribute over the Fermi 
surface in a more uniform way, is usually favored. 




FIG. 7: The spin configurations of the two-body states with 
a) J= 1 and J z = 0 and b) J = .J z = 0. The average magnetic 
dipolar interaction is attractive in a) but repulsive in b). From 
Ref. HJ. 


VI. THE J-TRIPLET COOPER PAIRING WITH 

THE MAGNETIC DIPOLAR INTERACTION 

The magnetic dipolar interaction also gives rise to 
novel Cooper pairing structures possessing the SO cou¬ 
pled nature^. Again, below, we use the simplest case of 
F = i as a prototype model to explore the exotic pair¬ 
ing structure of the magnetic dipolar interaction, which 
provides a novel and robust mechanism for the p-wave 
(L = 1) spin triplet (S = 1) Cooper pairing. It turns out 
that its pairing symmetry structure is markedly different 
from those in the celebrated p -wave 3 He A and B phases: 
the orbital angular momentum L and spin S' of a Cooper 
pair are coupled into the total angular momentum J = 1, 
and thus dubbed as the J-triplet pairing. In comparison, 
the 3 He B phase is isotropic in which J = 0; while, the 
A phase is anisotropic in which J is not well-defined^. 

Even within the J-triplet sector, there are still com¬ 
peting instabilities regarding to different possibilities of 
J z ’s: the helical polar state (J z = 0) which maintains TR 
symmetry, and the axial state (J z = ±1) which breaks 
TR symmetry. The helical polar state exhibits a pair 
of nodes around which the quasi-particle wavefunction 
exhibits winding numbers ±1, thus it is a gapless Dirac 
pairing. It is a TR invariant generalization of the 3 He-A 
phase with coupled spin and orbital degrees of freedom. 
This state was also proposed before, in the context of su- 
perfluid 3 He as an intermediate phase between the 3 He- 
B phase and the normal stated. In contrast, the axial 
pairing state exhibits a Weyl type node around which 
the winding number is 2, and its low energy spectrum is 
quadratic. 

In this section, we review the SO coupled pairing struc¬ 
ture of the magnetic dipolar fermions. An intuitive real 
space picture for the J-triplet pairing is presented in Sect. 
IVI Al The partial-wave analysis for the pairing eigenval¬ 
ues is performed in Sect. IVI Bl The structure of the nodal 
Bogoliubov quasi-particles is given in Sect. IVI Cl 


A. The real space picture for the J-triplet pairing 

We perform an intuitive real space analysis based on a 
simple two-body problem. It can be shown that the mag¬ 
netic dipolar interaction between two spin--f fermions 
vanishes in the total spin singlet channel, and it only 
exists in the total spin triplet channel. Naturally, the 
magnetic dipolar interaction leads to the triplet pair¬ 
ing. Let us introduce a characteristic length scale at 
which the kinetic energy equals the interaction energy 
at adp = mg 2 F g 2 B /h 2 where m is the fermion mass. Since 
our purpose is to seek the most attractive angular partial- 
wave channel, without loss of generality, we can fix the 
inter-particle distance at a,d P ■ In the spin triplet channel, 
the lowest value of orbital angular momentum is p-wave, 
and thus there are 3x3 = 9 states with L = S = 1. 
According to the total angular momentum J, they can 
be classified into J = 0,1 and 2. The interaction energies 
can be easily diagonalized in each sector as 

E 0 = Ed p , Ei = —~Ed p , E 2 = —Ed p , ( 96 ) 

respectively, where Ed p = g^pB/^dp- Only the channel 
with J = 1 can support bound states, which is also con¬ 
firmed from the momentum space analysis below in Sect. 

ETU 

The reason why the J-triplet channel is the most at¬ 
tractive one is because of its particular pairing spin con¬ 
figuration, which shows a dominant “head-to-tail” con¬ 
figuration and thus the interaction is attractive. More 
precisely, let us denote the spin wavefunctions and the 
relative orbital wavefunctions as \p and Pp(^)i which 
satisfy 

{e M • (Si + S 2 )} X p = 0, (e„ ■ L)pp(Q) = 0, (97) 

for p = x,y, and z. In other words, \p an d Pp(&) are 
polar eigenstates of total spin and the relative orbital 
angular momentum, respectively. In the sector of J = 1, 
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we define the SO coupled polar state satisfying 

(e M • J)<^ = 0, (98) 

which can be expressed as <fyj(fl) = -^£nvXXvPx{H)- I n 
particular, for the state with J z = 0, it can be further 
expressed as 


The interaction matrix element in Eq. 11021 is defined 


as 


Vs*S' z (k-,k') = l J2 {lS z \ttP){lS',\^Pr 


'2 2 ' 


'2 2 ' 


<y.pp' a.' 

X { VaP,P'ot' {k ^ ) V a p,p, a '(k H - & ) }• 


M&) 



sin0{|ai p )i 


\ote p )2 + |/3gp )l |/?e p )2 } , (99) 


The spherical harmonics decomposition of Vs z s’ z (k;k') 
can be formulated as 


where, e p = x cos </> + y sin <jx\ \a e ) and | j3 e ) are eigen¬ 
states of e p ■ a with eigenvalues of ±1, respectively. This 
shows the “head-to-tail” configuration in Fig. |T] (a), and 
thus, the corresponding interaction is attractive. 

In contrast, the eigenstate of J = 0 shows the “head- 
to-head” configuration as shown in Fig. [7] (6), and thus 
the interaction in such a state is repulsive. It is expressed 
as </> 0 (fl) = ^{|an)i|/3fi) 2 + |/3n)i|an) 2 }, where |an) 

and \/3q) are eigenstates of id • a with the eigenvalues ±1, 
respectively. 


B. The momentum space partial-wave analysis 


In this part, we review the partial-wave analysis in mo¬ 
mentum space, which arrives at the same pairing symme¬ 
try as that obtained through the real space analysis. 

After the mean-field decomposition, the pairing Hamil¬ 
tonian of the magnetic dipolar system becomes 


Hmf — 


d 3 k 

(2^ 


&(k) 


mi A a/3 ( k ) 

Afyfy/o -mi 


vfyfc), 

( 100 ) 


where ’F(fc) = (ipf(k), ^(k), i>\{~k), ^{(-fc)) T ; £(fc) = 
e(k) — [x. The Bogoliubov quasiparticle spectra become 


Ei(k) = yCfc + A?(fc) where i = 1,2 and Xf(k) are the 
eigenvalues of the positive-definite Hermitian matrix de¬ 
fined as Afyfc)A(fc). 

The pairing matrix A a p is defined as 


A a/3 = ^(15 s |iai/3)*A S2 

S z 


( 101 ) 


where ( lS z \^ah[i) is the Clebsch-Gordan coefficient for 
two spin-A states to form the spin triplet. A g z satisfies 
the mean-field gap function as 

_ i ,_ r J3U _ _ i 

i v ' 

x A s pk'), (102) 

where A',;(fc') = tanh[^A i (fc , )]/[2A,(/c')], and the integral 
in Eq. 11021 is already normalized following the standard 
procedure^. 


e 

LM,L'M' 


VLMS^L'M’S'YLM^k) 

x 1W(%)> (103) 


where L = L', or, L = I! ± 2, and L,L' are odd in¬ 
tegers. The expressions of the dimensionless matrix el¬ 
ements Vlms z \L'M'S' z are the same as those in Eq. [50] 
except an overall minus sign and a numeric factor. 

The free energy can be calculated as 




d 3 k , r_ . P E kAl 


In [2 cosh ■ 


S Z ,S‘ 

X As-(fc)}, 


(27r) 3 
d 3 k d 3 k' 

(2n) 3 (2n) 3 


{a imvsUm) 


(104) 


where V s g, (k;k') is the inverse of the interaction matrix 
defined as 


,_, /• ^3;./ _ _ _ 


In order to analyze the pairing eigenvalues, the gap 
equation is linearized around T c for states around the 
Fermi surface. The total angular momentum quantum 
number J is employed to classify the eigen-gap functions 
denoted as <p‘g JJz (k), in which the index a is used to 
distinguish different channels with the same value of J. 
The eigen-equation for <t>g’ JJz (k) is 


N, 


J V S*S' Z C k ; k')(/) a s ’/ Jz (k’) = w a j(j) a ^ Jz ( k ), (105) 


where Wj are dimensionless eigenvalues; k, k! are at the 
Fermi surface. Employing the spherical harmonics de¬ 
composition of Eq. 11031 the most negative eigenvalue is 
calculated lying the channel of J = L = 1 as w J=1 = 
—27rA m . All other negative eigenvalues are significantly 
smaller. Thus, the dominant pairing channel remains 
in the J-triplet channel in the weak coupling theory in 
agreement with the real space analysis. 

For later convenience, the pairing matrix in the J = 1 
sector is represented as 

A^(fc) = Y e ^x (fycr A - k\<j v ) al 3, (106) 
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COS0k 

FIG. 8: The angular distribution of the gap function |A(fc)| 2 
v.s. cos 9k in the helical polar pairing state (the red line) and 
the axial pairing state (the black lines). From Ref. [§oJ. 


for p = x,y and z. It represents a pairing symmetry 
whose angular momentum projection along the direction 
of is zero, i.e., it is an eigenstate of • J with the 
zero eigenvalue. 


C. Helical polar pairing and chiral axial pairing 

In the sector of J = 1, there are still two non-equivalent 
pairing possibilities: J z = 0, or J z = ±1, whose pairing 
matrices are A z a p(k) and ^j{A“ / 3 (A:)±iA^(/c)}, respec¬ 
tively. Based on the GL analysis up to the quartic or¬ 
der of the pairing order parameter, it can be proved that 
these two are the only non-equivalent pairing symmetries 
under 3D rotations. Right at T c , the Ginzburg-Landau 
free energy can be linearized, and these two instabilities 
are degenerate, while this degeneracy is lifted below T c 
due to the non-linearity of the Ginzburg-Landau free en¬ 
ergy. 

In quantum mechanics, if a system possesses rotation 
symmetry, of course, the eigenstates of its energy in the 
J = 1 sector and all of their superpositions are degener¬ 
ate because of the linearity of quantum mechanics. How¬ 
ever, the index J labeling the pairing order parameter is 
not the angular momentum of the many-body eigenstate. 
The description of a many-body system in terms of or¬ 
der parameters is a great simplification by only keeping 
a very limit but essential amount of degrees of freedom. 
A price to pay is that the description in terms of or¬ 
der parameters, say, the Landau-Ginzburg free energy, is 
non-linear even though quantum mechanics remains lin¬ 
ear. This is the reason why the J z = 0 and J z = ±1 
in principle are non-equivalent and the superposition law 
does not hold for Cooper pairing symmetries. The two 
pairing patterns with J z = 1 and J z = — 1 are equivalent 
to each other which can be connected either by a rotation 
or by the TR transformation. 


1. Helical polar pairing 

The pairing symmetry with J z = 0 is also called the 
J-polar pairing, characterized by the following pairing 


matrix A^ = \\A\[k y ai - k x a 2 )i<J2]ai3 as 

k x — iky 0 

0 k x + iky J ’ 



(107) 


which is a unitary pairing and preserves TR symmetry. 
The spin-up fermions are paired in the p x — ip y symme¬ 
try, while, the spin-down fermions are paired with the 
Px T iPy symmetry. Thus, it is equivalent to a helical ver¬ 
sion of the 3 He-A phase. In the 3 He-A phase, the pairing 
symmetry defines an orbital angular momentum direc¬ 
tion represented by the Z-vector. Here, the z-axis plays a 
similar role, but due to the TR symmetry, it represents 
the bi-direction of the polar axis, and thus it is no long a 
vector but a director. Thus if we rotation the polar pair¬ 
ing around any axis in the a,’y-plane for 180°, then the 
system returns to itself up to a global phase difference. 

The Bogoliubov quasiparticle spectra are degener¬ 
ate for two different spin configurations as E? = 

+ \AP l (k)\ 2 with the anisotropic gap function 

\A pl (k)\ 2 = j|A| 2 sin 2 0fc depicted in Fig. [5] The effec¬ 
tive linearized Hamiltonian around the north and south 
poles can be combined as 


H n =ipi(k) 


Vf{k z - kf) k x - iky 
k x + ik y -Vf(k z -kf) 


(108) 


where i/> T (k) = (^-j-(fc), ip\(—k))- Eq. 11081 describes a 
Weyl fermion with a unit winding number. Combining 
the effective Hamiltonian for the spin down sector which 
is also a Weyl fermion but with an opposite chirality, the 
low energy Bogoliubov spectrum is 3D gapless Dirac like. 
If the quantization axis for the polar pairing is rotated 
away from the z-axis, then the spin quantization axis for 
the Bogoliubov quasiparticles should also be transformed 
accordingly. 


2. The chiral axial pairing 


Without loss of generality, we pick up the J z = 1 pair¬ 
ing, and the result of the state of J, = — 1 can be obtained 
by performing TR transformation. The pairing matrix, 
A a /3 = ^\^\{k z {cri+ia 2 )i<T 2 +ia 3 a 2 (k x +iky)} a p, takes 
the form of 


A %$) = ^|A| 


L 2 


2 ( n 'I “T bFVy ) 

2^k x T iky') 0 


2 i,k X + iky') 


, (109) 


thus, this is a non-unitary pairing state. The Bogoliubov 
quasiparticle spectra have two non-degenerate branches 
with anisotropic dispersion relations as E^ x {k) = 

+ |Af x (fc)| 2 , with 


|A“ x (fe )| 2 = ^|A | 2 (1 ± cos #*,) 2 ( 110 ) 


for i = 1 and 2, respectively, as depicted in Fig. [8] 
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Considering the coupling between k and —fc, we can 
combine the north and south poles together into a four 
component spinor. The above energy spectra show that 
two of them are gapped, while, another two of them are 
gapless forming a two-component Weyl spinor. Different 
from the usual linear dispersion of Weyl fermions, the dis¬ 
persions are quadratic with respect to the the transverse 



A natural question is: which pairing is energetically 
more stable? For the case of 3 1 He, at the mean-field level, 
the H-phase is always more stable than the A-phase^ 
because the gap function of the B phase is uniform over 
the Fermi surface. This can be intuitively understood as 
follows: the BCS mean-field free energy density in mo¬ 
mentum space F\\ A(fc)| 2 3 4 ] can be viewed as a functional 
of the gap function |A(fc)| 2 . Let us consider the pairing 
on the Fermi surface with the constraint of / dflfc|A(A;)| 2 
fixed as a constant, and then to minimize the free en¬ 
ergy. Usually, the non-linearity of the free energy favors 
a distribution of |A*,| 2 as uniform as possible. Here, the 
situation is quite subtle. As shown in Fig. 0 for both 
cases, the pairing gap function distributions on the Fermi 
surface are non-uniform. Nevertheless, the distribution 
of the unitary polar pairing is more uniform than that 
of the axial pairing. Naturally, we expect that, at the 
mean-field level, the polar pairing wins. This has been 
numerically confirmed by comparing the BCS mean-field 
free energies of both pairings in Ref. [89| . However, we 
need to bear in mind that this conclusion is only valid at 
the mean-field level. We cannot rule out the possibility 
that certain strong coupling effects can stabilize the axial 
state. In fact, even in the 3 He case, the A phase wins in 
certain parameter regime in which the strong correlation 
spin feedback effect dominate s 56 ! 88 . 


optical lattice structure. For particles in the same lattice 
site, the inter-particle distance is at the order of lOOwn, 
which can enhance the interaction energy scale to 0.6/zA'. 
Even the s-orbital band can hold up to 2 F + 1 fermions 
per site. We expect that this system can exhibit a variety 
of even more exotic many-body physics to be explored. 


VII. CONCLUSIONS 

We have briefly reviewed the novel many-body physics 
of both electric and magnetic dipolar fermion systems 
focusing on the aspect of unconventional symmetries. 
The electric dipolar interaction is characterized by its 
<j r 2 _ 32 2 -anisotropy, which leads to anisotropic Fermi liq¬ 
uid properties in the particle-hole channel and a robust 
mechanism of spin-triplet p z -wave Cooper pairing in the 
particle-particle channel. The competition and coex¬ 
istence between the singlet and triplet Cooper pairing 
leads to a novel mechanism of TR symmetry breaking 
pairing. The magnetic dipolar interaction manifests its 
SO coupled nature in unpolarized fermion systems. Its 
Fermi liquid theory is SO coupled whose collective zero 
sound mode exhibits an oscillation of a topological non¬ 
trivial spin structure over the Fermi surface. The mag¬ 
netic dipolar interaction also lead to a SO coupled p -wave 
spin triplet Cooper pairing state with the total spin of a 
Cooper pair J = 1. This is a novel pairing symmetry 
different from that in both 3 He-H and A phases, whose 
Bogoliubov quasiparticles exhibit nodal helical Dirac or 
chiral Weyl spectra. 
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